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1. Introduction and motivation 

Chiral perturbation theory (#PT) ^ and lattice QCD are powerful tools to extract quantities at 
low energy which are relevant for QCD. Lattice QCD simulations can provide the determination of 
the Low Energy Constants (LECs) of #PT from first principles calculations. A particular situation 
is found when one enters the £ -regime [0], where 
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(1.1) 



There are analytical formulae from #PT that describe the behavior of physical quantities in the 
£ -regime as a function of the volume and the quark mass. These formulae are parameterized by the 
infinite volume LECs of the effective chiral Lagrangian. Therefore we can extract physically rele- 
vant information even from the unphysical £-regime. As a peculiarity of the £-regime, since one is 
working in small boxes (L> 1.1 fm- • • 1.5 fm), observables depend significantly on the topological 
sector, and predictions exist for expectation values in specific sectors [Qj. For the parameters that 
have been used in the £-regime simulations, it would be of particular interest to collect large sets 
of configurations 1 with an index | v | > 0. 

In general, however, it is not obvious to define topological sectors on the lattice. A neat 
definition exists with overlap fermions satisfying the Ginsparg-Wilson relation [pi, 0] 
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where ju > 1 is a mass parameter in the operator Q. Since D m has exact zero modes with definite 
chiralities [|6[ [7|, one can use the index as a definition of the topological charge due to the 
Atiyah-Singer theorem. 

While overlap fermions define the fermion sector of lattice QCD, if one insists on exact lattice 
chiral symmetry, the lattice gauge action of QCD can be constructed in many ways One 
requires the naive continuum limit for all lattice gauge actions to coincide, in which case they fall 
into the same universality class. The simplest formulation is the Wilson plaquette action 



S W [U] = PY,MUp) , S P (U P ) = 1 - -Re TtU p , 



(1.3) 



where the sum over P runs over all plaquettes. 



For the overlap operator (1.2), the topological transitions are excluded under continuous defor- 



mations if all the plaquette variables Up in the configurations involved obey the inequality [10, 11] 

2 J_ 

5d{d-\) 30 

1 



Sp(U P 

£ 



< £ 



(for all P). Later on H. Neuberger showed a more tolerant bound [ |12| ] 

(i+i/V2)d(d i) — 2S3 ■ F° r tn * s admissibility condition, the exponential locality of the Ginsparg- 
Wilson fermions is guaranteed [[TT|], and the topological charge is conserved rigorously. However, 
it is very difficult to realize this condition in a numerical simulation. For practical purposes we 
have to relax £ to larger values than the analytical bounds. 



'The topologically neutral sector is problematic due to the frequent appearance of very small Dirac eigenvalues, 
which leads to strong spikes in the Monte Carlo histories of correlation functions Mm. 
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2. Proposal of the gauge actions 

We now describe a number of non-standard lattice gauge actions, which suppress the unwanted 
small plaquette values leading to transitions between different topological sectors. The naive con- 



tinuum limit of these actions coincides with that of Wilson action (1.3). Based on the proposal by 



M. Luscher [10], we consider topology conserving gauge actions [14] 



[l-S P (Up)/e]» 



for Sp(Up) < e and +oo otherwise, 
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fi$(U P ) = S p (Up) ■ exp{Sp(U P )"/e} , 



Sp(Up) + -Sp(U P ) n , 
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In simulations, the action ( |2.ip with n = 1 was first used in the Schwinger model by Fukaya 
and Onogi [13]. They set £ = 1, i.e. far above the theoretical value of about 0.29, but they still 



observed topological stability over hundreds of configurations. Here we investigate this type and 
its extension to the gauge action of lattice QCD. 

We show the basic properties of these actions. The left plot of fig. |T] is the histogram of 
plaquette values for the action S^ =l , compared to the Wilson action, on a 4 4 lattice. In the topology 
conserving gauge action, the occurrence of very small plaquette values is drastically suppressed. 
The right plot of fig. [I] shows the ratio of the force in the HMC algorithm between the Wilson 
plaquette action Sp and topology conserving gauge actions. The action Sff makes a sharp wall 
continuous, but shows the same behavior as the Wilson action over a wide range. 





Figure 1: Left: Histograms of Sp for the topology conserving gauge action, compared to the Wilson action, 
on a 4 4 lattice. Right: Ratio of the HMC force between the Wilson action Sp and modified actions. 

For the generation of the configurations we used a local HMC algorithm. Since S^ is non- 
linear for the link variables, the heatbath and over-relaxation algorithms are not straightforwardly 
applicable. We use a 16 4 lattice and measure physical quantities every 50 trajectories. See table |] 
about the acceptance rate of local HMC. At 1/e = 1.64, the acceptance rate is about 65% at the 
Molecular dynamics step d% = 0.1. However the acceptance improves as dx decreases. 



3. Monte Carlo history of Qtop' Stability of the topological charge. 

In fig. ^| we show the Monte Carlo history of the topological charge Q top , which is evaluated 
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act.=Wilson, p=6.18, 1/e=0.00, dt=0.1 




act:Hyp,p=0.8, 1/e=1.25, Q , =1, dt=0.1 
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Figure 2: MC history of Q top for S w (left) and S^? ag „ =1 (right). 



by the cooling method 2 , for the Wilson action and a modified gauge action. The topological charge 
in 1/e <C 20.5 is not always conserved. However, the changes of Q top are suppressed. Here we 
define the quantity representing the stability of the topological charge: f top =(the number of jumps 
of Qtop) I (the number of trajectories), see also table [IJ The value f top for the modified gauge actions 
is 10 times smaller than for the Wilson action and thus these actions clearly stabilize the topology 
charge. Also the autocorrelation of the plaquette x plaq is getting shorter as 1 /e increases. 



4. The static potential, the physical scale and lattice artifacts 

In this section we evaluate the physical scale and the lattice artifacts for various parameters of 
/3 and l/e in the modified gauge actions. The most established method of setting a scale in pure 
gauge theory is the measurement of the static potential and the force at relatively large distances, 



extracted from Wilson loops, see e.g. Ref. [15]. 



Results are shown in the left plot of fig. ||. A way to check the lattice artifacts [15] is to 
compare the short distance force at finite lattice spacing with the one extrapolated to the continuum 
limit r 2 F(r/r )\ c 

4(r/ro) = ^W^«. (4,, 
r z F[r/ro)\ c 

As seen in this figure, for increasing l/e, the discretization errors increase as well. Typically the 
artifacts of S^ =l at l/e £1 are about 10%. Even in l/e = 1.64, the lattice artifact is less than 
15%. This is comparable with the Iwasaki and DBW2 gauge actions. Note that the static potential 
at distance r ~ r$ has not been determined with high precision (this is reflected in large uncertainties 
on the quantities A(r/rn)). 

5. Condition number for the overlap operator 

We also evaluate the effect of topology conserving gauge actions on the condition number of 



2 The cooling charges agree in practically all cases with the overlap indices in the range jl = 1.3, ...,1.6, as we 
verified for a subset of the configurations. 
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Figure 3: Left: Lattice artifact by the standard estimate for the action S^ =l . Right: The overlap condition 
numbers c^.-.c^x at jj. = 1.6. 



the overlap operator (1.2). We observe the condition number defined as C„ = -f^ where X max is 
the maximum eigenvalue and A„ is the rc'th eigenvalue of the kernel Q 2 in the overlap operator. 
C„ is relevant if n — 1 modes of Q 2 are projected out. Fig. ^ (right plot) shows the behavior of 
the condition numbers C2, ■ ■ • ,C2\ at jU = 1.6 and about r^/a ~ 7.0. We see an improvement which 
increases if only a few modes are projected out and if 1/e grows. Hence the computation of 
Q/ \[&~ in the overlap operator gets fast by using topology conserving gauge actions. 



6. Summary 

Simulations in the £-regime of chiral perturbation theory using the overlap operator can be 
simplified when appropriate choices of the gauge actions are made. We tested a number of gauge 
actions that suppress small plaquette values and investigated the properties of these actions fll6|]. 
Our results can be summarized as follows, see also table [[]. 



The actions of eqs. ( gjj , [2.2| , |2.3[ ) stabilize the topological charge and could therefore be 
profitable in QCD simulations in the £-regime. 



• The actions of eqs. (^2J, [2.3D are conceptually clean, have a positive transfer matrix and show 
lattice artifacts that are acceptable and compare well with those of the Iwasaki and DBW2 
gauge actions. 

• The condition number of the kernel of the overlap operator is improved by using the topology 
conserving gauge actions when compared to the standard Wilson plaquette action. This 
speeds up the simulations of overlap fermions in the quenched approximation and will also 
help in dynamical simulations using the (global) HMC algorithm. 
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P 


ro/a 


Pw 


dr 


T plaq 


/top 


acc. rate 


0. 


6.19 


7.14(3) 


6.19 


0.1 


7(1) 


2.2e-2 


> 99 % 


1. 


1.5 


6.6(2) 


6.13(2) 


0.1 


2.2(1) 


2.4e-3 


> 99 % 


1. 


1.5 


6.6(2) 


6.13(2) 


0.01 


2.2(1) 


3.2e-3 


> 99 % 


1.18 


1. 


7.2(2) 


6.18(2) 


0.1 


1.2(1) 


1.6e-3 


> 99 % 


1.18 


1. 


7.2(2) 


6.18(2) 


0.02/0.01 


1.3(1) 


1.4e-3 


> 99 % 


1.25 


0.8 


7.0(1) 


6.17(1) 


0.1 


1.1(1) 


2.5e-3 


> 99 % 


1.52 


0.3 


7.3(4) 


6.19(4) 


0.1 


0.8(1) 


9.4e-4 


«95 % 


1.64 


0.1 


6.8(3) 


6.15(3) 


0.1 


1.0(1) 


7.0e-4 


ss 65 % 


1.64 


0.1 






0.05 


0.7(1) 


2.3e-3 


«78 % 


1.64 


0.1 






0.025 


0.6(1) 


3.5e-3 


«93 % 



Table 1: Results for Sff at various values of e and j3, on a 16 4 lattice, fiw is the j3 value of Wilson action 
corresponding to the measured ratio ro/a. x plaq is the autocorrelation time of the plaquette value. 
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